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ABSTRACT 


Shock  tube  test  time  limitation  due  to  the  premature  arrival  of  the 

contact  surface  is  analytically  investigated  for  wholly  turbulent 

wall  boundary  layers.  The  results  are  compared  with  those  for 

wholly  laminar  wall  boundary  layers.  It  is  found  that  for  a  given 

shock  Mach  number,  M  ,  the  maximum  possible  test  time  (in  a 

®  5/4  1/4  2 

long  shock  tube)  varies  as  d  p  and  d  p  for  the  turbulent 

00  00 

and  laminar  cases,  respectively  (d  =  tube  diam,  p^  ==  initial 

pressure).  For  3  <  <  8  in  air  or  argon,  it  is  found  that  the 

turbulent  boundary  layer  theory  for  maximum  test  time  applies 

for  dp^  ^  4  to  10  (air)  and  ^  2  to  10  (argon),  where  d  is  in 

inches,  p^  is  in  cm  Hg.  Similarly,  for  3  <  M^  <  8,  the  laminar 

theory  applies  for  dp^  ^0*3  (air)  and  dp  0*  6  to  1  (argon). 

When  %  ^9  turbulent  theory  for  both  air  and  argon  indicates 

test  times  of  about  1/2  to  1/4  the  ideal  vadue  for  =  43  to  150» 

respectively  (x^  =  length  of  low  pressure  section).  Higher  values 

of  dp  result  in  more  test  time.  When  dp  ^  0.  5,  laminar  theory 
00  00 

indicates  about  1/2  ideal  test  time  for  x  /d^  100.  Lower  dp _ 

8  *^00 

reduces  test  time.  Working  curves  are  presented  for  more  accu¬ 
rate  estimates  of  test  time  in  specific  cases. 
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I.  INTRODUCTION 


In  an  ideal  invitcid  shock  tube,  the  separation  distance  between  shock 
and  contact  surface  increases  linearly  with  distance  from  the  diaphragm 
(Fig.  1).  In  a  real  shock  tube,  however,  the  wall  boundary  layer  between  the 
shock  and  contact  surface  acts  as  an  aerodynamic  sink  and  absorbs  mass 
from  this  region.  This  causes  the  contact  surface  to  accelerate  and  the 
shock  to  decelerate  and  reduces  the  separation  distance,  i,  below  the  ideal 
value.  The  separation  distance  approaches  a  limiting  (maximum)  value, 

At  this  limiting  condition,  the  shock  and  contact  surface  both  move  with  the 
same  velocity  (Fig.  1).  The  limiting  separation  distance  is  that  at  which  the 
mass  flow  entering  the  shock  equals  the  boundary  layer  mass  flow  moving 
past  the  contact  surface. 

Separation  distance  imposes  an  upper  bound  on  the  test  time  in  shock 
tubes.  (Nonuniformity  of  the  flow  between  shock  and  contact  surface  may 
further  reduce  test  time.  The  amount  of  flow  nonuniformity  that  can  be 
tolerated  depends  on  the  nature  of  the  experiment  and  the  instrumentation.) 

It  is  thus  important  to  know  separation  distance  as  a  function  of  distance  from 
the  diaphragm  in  order  to  estimate  test  times  in  shock  tubes. 

The  effect  of  a  laminar  wall  boundary  layer  on  separation  distance  has 
been  studied  experimentally  by  Duff  (Ref.  1)  and  both  experimentally  and 
analytically  by  Roshko  (Ref.  2)  and  Hooker  (Ref.  3).  Some  aspects  of  the 
analyses  presented  in  Refs.  2  and  3  are  modified  in  Ref.  4.  These  references 
can  be  used  to  estimate  test  time  in  low  pressure  shock  tubes  (roughly,  those 
shock  tubes  with  initial  pressures  of  the  order  of  1  mm  Hg  or  less,  currently 
being  used  to  study  dissociation,  ionisation,  and  other  rate  phenomena).  The 
test  time  limitation  in  low  pressure  shock  tubes  has  received  considerable 
attention  because  the  upper  limit  on  test  time  is  proportional  to  d^p 
(Ref.  2)  and  very  low  d  p^cen  result  in  virtually  no  usable  test  time 
(d  s  tube  diam,  ■  Initial  pressure). 
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Fig.  1.  Boundary  Layer  Effect  on  ^ock  Tube  Flow 


The  effect  of  a  turbulent  wall  boundary  layer  on  separation  distance  has 
been  treated  only  by  Anderson  (Ref.  5).  It  is  interesting  to  note  that  Anderson 
was  the  first  to  blame  less  than  ideal  shock  tube  test  time  on  the  premature 
arrival  of  the  contact  surface.  Anderson  showed,  for  two  practical  examples, 
that  mass  loss  to  a  turbulent  wall  boundary  layer  reduced  the  theoretical  test 
time  by  about  a  half.  This  was  in  agreement  with  a  rule-of-thumb  presented 
in  Ref.  6  for  shock  tubes  used  to  study  re-entry  heat  transfer.  In  Ref.  6,  it 
was  stated  that  40  to  70  percent  of  theoretical  test  time  was  generally  ob¬ 
tained,  and  BO  percent  was  proposed  as  a  mean  value.  (The  test  time 
referred  to  was  the  period  of  essentially  uniform  conditions.  The  reduction 
from  the  ideal  value  was  attributed,  in  Ref.  6,  to  mixing  at  the  interface  and 
flow  nonuniformity  associated  with  shock  attenuation  rather  than  to  the  pre¬ 
mature  arrival  of  the  contact  surface).  Since  test  time  did  not  appear  to 
represent  a  serious  problem  in  re-entry  heat  transfer  shock  tubes  studies, 
the  work  of  Anderson  did  not  find  widespread  application. 

As  previously  noted,  the  effect  of  a  laminar  boundary  layer  on  shock 
tube  test  time  is  now  receiving  considerable  attention  (Refs.  1-4).  It  is 
therefore  felt  that  it  is  worthwhile  to  develop  the  turbulent  boundary  layer 
case  beyond  the  preliminary  work  of  Anderson.  Anderson  presented  a  first 
estimate  of  the  variation  of  separation  distance  with  distance  from  the 
diaphragm  for  a  particular  initial  pressure  (10  cm  Hg)  and  two  tube  diameters 
(HI/ 2  and  4  in.).  In  the  present  paper,  the  problem  is  formulated  more 
accurately,  and  the  results  are  presented  in  nondimensicnal  form  so  as  to  be 
applicable  for  arbitrary  initial  pressures  and  tube  diameters.  Criteria  are 
established  to  define  when  turbulent  wall  boundary  layer  theory  is  applicable 
and  when  laminar  wall  boundary  layer  theory  is  applicable.  The  development 
here  is  parallel  to  that  used  in  Ref.  4  for  the  laminar  case. 

The  present  solution  is  primarily  intended  for  the  case  where  the  wall 
boundary  layer  introduces  sizable  reductions  in  test  time.  When  the  wall 
boundary  layer  introduces  only  small  perturbations  of  the  ideal  test  time,  the 
linearized  methods  of  Refs.  15  and  16  are  more  appropriate. 
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(b)  Uniform  Free  Stream 


Fig.  2.  Flow  Between  Shock  and  Contact  Surface  in 
Shock>Stationary  Coordinate  Syatem 
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n.  MAXIMUM  SEPARATION  DISTANCE 


Due  to  the  wall  boundary  layer,  the  separation  distance  between  shock 
and  contact  surface  approaches  a  maximum  value  such  that  the  mass  flow 
through  the  shock  eqixals  the  mass  flow  moving  past  the  contact  surface. 

’‘ihin  phenomenon  has  been  observed  for  laminar  wall  boundary  layers  (Refs. 
1^3)  and  may  also  be  assumed  to  occur,  for  suificiently  long  shock  tubes,  in 
the  case  of  a  turbulent  wall  boundary  layer.  When  the  maximum  separation 
distance  is  reached,  both  shock  and  contact  surface  move  with  the  same 
velocity.  The  flow  between  the  shock  and  contact  surface  may  then  be 
viewed  as  steady  in  a  coordinate  system  in  which  the  shock  is  stationary.  In 
this  shock- stationary  coordinate  system,  the  wall  moves  with  velocity  u^ 
(which  equals  the  shock  velocity  Ug  the  laboratory  system).  This  steady 
flow  is  investigated  in  this  section  with  the  primary  object  of  determining  the 
maximum  separation  distance  between  the  shock  and  the  contact  surface.  The 
variation  of  separation  distance  with  distance  from  the  diaphragm  is  treated 
in  Section  III. 

The  steady  flow  is  illustrated  in  Fig.  2a.  The  shock  is  located  at 
1  =  0,  and  the  free  stream  portion  of  the  contact  surface  at  i  =  The  flow 

upstream  of  the  shock  is  denoted  by  subscript  oo  and  moves  with  velocity  u^, 
as  does  the  wall.  Free  stream  conditions  between  the  shock  and  the  contact 
surface  are  denoted  by  subscript  e.  Free  stream  conditions  directly  down¬ 
stream  of  the  shock  have  the  additional  subscript  o.  The  percentage  of  mass 
^ow  in  the  boundary  layer  increases  with  i  such  that  all  the  mass  flow  is  in 
the  boundary  layer  at  and  the  free  stream  is  stationary  at  that  location. 


e 

For  laminar  boundary  layers,  the  shock  and  contact  surface  velocities  tend 
to  remain  constant  after  the  maximum  separation  distance  is  reached 
(Ref.  1).  For  turbulent  wall  boundary  layers,  the  shock  and  contact  surface 
may  both  continue  to  decelerate  after  the  maximum  separation  is  reached. 
This  deceleration  depends  on  the  driver  gas.  If  the  shock  continues  to 
decelerate,  then  the  present  steady  solution  is  a  quasi -steady  solution. 


.5. 


The  physical  parameters  upon  which  depends  can  be  found  as 
follows.  The  flow  rate  through  the  shock,  th^,  is 


where  A  is  the  cross-sectional  area  of  the  tube.  If  it  is  assumed  that  the 
boundary  layer  is  thin  relative  to  the  tube  radius,  the  mass  flow  in  the 
boundary  layer  at  the  contact  surface  location,  m^,  can  be  characterized  by 


m 

c 


Lp  (u 
'^w,  o  w 


u 

e, 


)6 


(lb) 


where 


(Ic) 


Here,  L  is  the  perimeter  of  the  tube;  is  a  characteristic  boundary  layer 

displacement  thickness  at  i  ;  p  and  u  -  u  are  characteristic  densities 

m  ^w,  o  w  e,  o 

and  velocities,  respectively;  p  is  a  constant;  and  n  =  1/2  for  laminar  and  1/5 
for  turbulent  boundary  layers.  Equating  Eqs.  (la)  and  (lb)  yields 


i 


1  -n 
m 


e,  o 


4P  p 


w,  o 


u  /u 

e,  o  /  w 

u  -  u  V 
w  e,  o\ 


-  u 
_ ^ 


V 

w,  o 


(2) 


where  d  =  4A/L  is  the  hydraulic  diameter  of  the  tube. 

Assume  that  the  temperature  upstream  of  the  shock  is  at  a  standard 
condition  so  that  T_  =  T  ,,  a  =  a  and  u  =  u  . .  Also,  assume  that  the 

OO  St  00  St  00  St 
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Equation  (2)  can 


wall  remains  at  its  initial  temperature  so  that 
now  be  put  in  the  form 


=  T 


St’ 


(3) 


where  W  =  u  /u  -  p  /p  ;  M  =  u  /a  ;  and  p  ^  is  a  standard  pressure, 
w'  e,o  ^e,o  ^oo  s  w'  oo  ^st  ^ 

usually  taken  to  be  one  atmosphere.  For  n  -  1/2,  these  equations  reduce  to 
those  presented  in  Ref.  4.  For  n=  1/5,  Eq,  (3)  becomes 


(4) 


The  right-hand  side  of  Eq.  (3)  depends  primarily  on  Mg.  For  a  given  Mg, 
i  ~  d^p  for  a  laminar  boundary  layer,  and  i  ~  d^/^  p^/^  for  a  turbu- 
lent  boundary  layer.  Therefore,  for  turbulent  boundary  layers,  is  less 
sensitive  to  variations  of  d  and  p^  than  is  the  case  for  laminar  boundary 
layers. 


Equation  (4)  does  not  yield  numerical  results  for  unless  an  accurate 
estimate  of  p  is  available.  Actually,  P  must  be  found  from  rather  than 
vice  versa.  That  is,  Eq.  (2)  is  taken  as  defining  P,  that  is, 


\ 


P  = 


d 


4i 


1-n 

m 


(5) 


and  is  found  as  accurately  as  possible  from  a  consideration  of  the 
boundary  layer  development  in  Fig.  2a.  A  first  estimate  for  P,  termed  P^, 
is  made  below  by  assuming  the  boundary  layer  to  develop  in  a  uniform 
external  stream.  An  improved  estimate,  termed  p^,  is  then  made  by 
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considering  a  variable  external  stream  and  employing  the  concept  of  local 
boundary  layer  similarity. 

A.  UNIFORM  EXTERNAL  FREE  STREAM  APPROXIMATION 

Boundary  layer  development  for  the  case  of  an  external  free  stream 
that  does  not  vary  with  i  is  illustrated  in  Fig.  2b  and  is  discussed  in 
Appendix  A.  Let  correspond  to  the  value  of  i  at  which  the  excess  mass 
flow  in  the  boundary  layer  equals  the  mass  flow  entering  through  the  shock. 
The  mass  flow  past  the  contact  surface  is  then 


m 


L(p  u  ) 

^  e  e  o 


m 


K  V-£. 

o  m  \a  -  u 
\  w  e 


(6) 


where  6*  is  the  displacement  thickness  at  f  and  K  is  a  function  defined  in 
m  mo 

Appendix  A.  Equating  Eqs.  (lb),  (Ic),  and  (6)  then  yields 


P 


o 


Pqo  W  ^ 

p  W  -  1  O 
e,  o 


(7) 


Values  of  P^  have  been  computed  for  y  =  7/ 5,  o*  =  0.  72  and  y  =  5/3,  <r  ~  0.  67 

using  the  turbulent  boundary  layer  theory  of  Appendix  A  together  with  the 

Sutherland  viscosity  law.  (y  -  ratio  of  specific  heats,  o'  -  Prandtl  number). 

Values  of  P^  have  also  been  computed  for  the  flow  behind  strong  shocks  in 

air  at  initial  pressures  of  0.  5  and  10  cm  Hg.  These  four  cases  are  further 

identified  in  Table  1  and  the  numerical  results  for  P  are  listed  in  Table  2. 

o 

These  values  of  P^  can  be  used  to  find  turbulent  boundary  layer  displacement 
thickness  (Eqs.  6,  7)  and  will  be  used  later  to  find  P^. 

It  is  expected  that  these  values  of  P^  will  overestimate  I particularly 
for  shock  Mach  numbers  that  are  not  large,  due  to  the  fact  that  the  relative 
velocity  between  the  wall  and  the  free  stream  increases  from  u^  -  u^  ^  at 
i  =  0  to  u  at  i  (compare  Figs.  2a  and  2b).  Hence* the  excess  mass  flow 

Vr  m 
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in  the  boundary  layer  will  be  greater  at  a  given  i  than  the  excess  obtained 

from  the  above  model.  This  will  result  in  smaller  i  and  larger  P  than 

m 

obtained  from  Eq  (7),  However,  for  very  strong  shocks,  where  u  is 

e ,  o 

small  relative  to  u^,  Eq.  (7)  should  give  accurate  results. 

B.  LOCAL  SIMILARITY  APPROXIMATION 

In  this  subsection,  the  streamwise  variation  of  free  stream  properties 
due  to  the  increase  in  boundary  layer  mass  flow  with  i  is  taken  into  account. 
The  development  of  the  boundary  layer  and  the  variation  in  free  stream 
properties  are  treated  simultaneously.  The  boundary  layer  growth  is  found 
by  assuming  that  at  each  station  it  is  similar  to  a  corresponding  boundary 
layer  developing  in  a  uniform  free  stream  behind  a  shock  moving  with  uniform 
velocity  (i.e.,  local  similarity). 

Since  the  flow  is  steady  (Fig.  2a),  the  net  mass  flow  through  the  shock 
equals  the  net  mass  flow  at  any  station  i  .  Thus 


A(p  u  )  =  Ap  u  + 

^e  e  o  ^e  e 


/<xy 

i.r(pu  - 


Pe'^e^  dy 


(8) 


In  Eq.  (8),  it  is  assumed  that  the  boundary  layer  thickness  is  small  compared 
with  d;  thus  the  integrand  is  nonzero  only  in  the  region  close  to  the  wall. 
Otherwise,  for  circular  tubes,  the  coefficient  L  would  have  to  be  replaced  by 
a  factor  ir(d  -  2y)  in  the  integrand  of  Eq.  (8).  Define 


5 


A  P  ^ 

4  ^e  e 

d  (p  u  ) 
^e  e  o 


(9) 


where  6^  is  the  boundary  layer  displacement  thickness  based  on  the  local 
free  stream 

6*  =  f°°{l  dy  (10) 

•'o  Pe  e 
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JL 


Fig.  i.  Boundary  Layer  Development  With  Uniform  Free 
Stream,  Wall  Velocity  and  Origin  at  i. 


Note  that  5  is  the  ratio  of  the  excess  mass  flow  through  the  boundary  layer  at 
i  to  the  mass  flow  through  the  shock.  Thus  6  varies  from  0  at  i  =  0  to  1  at 
i  =  Equations  (8)-(10)  then  give 


6  =: 


Pe% 

(Firj 

'^e  e  o 


(11) 


which  relates  the  free  stream  conditions  to  the  local  boundary  layer  displace¬ 
ment  thickness. 

The  concept  of  local  similarity  is  now  introduced.  It  is  assumed  that 
the  boundary  layer  profile  at  each  i  corresponds  to  the  profile  for  a  boundary 
layer  associated  with  a  uniform  free  stream  (equal  to  the  local  free  stream) 
and  a  wall  velocity  u^.  The  origin  of  this  fictitious  boundary  layer  is  at  i., 
which  is  initially  an  unknown  function  of  i  .  (See  Fig.  3.)  The  origin  i.  is 
chosen  such  that  the  excess  flow  the  boundary  layer  at  each  i  has  the 
correct  local  value.  sufficiently  large  that  variations  in  and  p^ 

can  be  ignored,  Eqs.  (11)  and  (A -11a)  become 

6  =  1  -  V  (12a) 

(12b) 

where  V  =  u  /u  and  H  is  a  function  of  the  local  free  stream  as  defined  in 
e  e,  o  e 

Eq.  (A-llb).  The  problem  is  now  to  solve  these  equations  simultaneously  to 
find  6  as  a  function  of  i.  The  value  of  i  at  6  =  1  defines  i  and  thus  pro- 
vides  an  improved  estimate  for  p. 


I 
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Assume  that  the  local  rate  of  growth  of  8  follows  the  local  similarity 

law  6  =  H  (i  -  i .)  *  (withH  and  i.  considered  constant  at  the  local  values), 
e  i  e  1 

This  is  the  approach  used  in  Ref.  4.  It  follows  that 


A5 

S 


^  (1-n) 


Ai 


f  -i. 


l/(l-n) 


=  (l-n)  ^ 


Af 


(13) 


For  n  =  1/2,  this  expression  is  the  same  as  that  derived  in  Ref.  4.  For 
n  =  1/5,  Eqs.  (13)  and  (12a)  give,  in  integral  form. 


_  5  r  *  (1  - 

'4Jv  „5/4 


dV 


(14) 


Substituting  Eq.  (A-llb)  for  and  considering  the  limit  V  =  0,  i  = 


yields 


(1  -V) 


1/4  W  +  BV/W-  l\^ 
W  +  B  Vw  -  V/ 


dV  = 


F 


(15) 


The  constant  B  is  defined  in  Appendix  A.  Equation  (15)  can  not  be  integrated 
in  closed  form.  For  W  >  2,  F  can  be  expressed  as 


F  = 


iB  i  ^  1 
B  +  W  W  -  1 

(5  +  4i)(l  -  W)^"^ 


(16) 
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so  that  F  -►  1,  from  below,  as  W  -♦oo.  Equation  (15)  has  been  numerically 
integrated  for  B  =  7/3  and  2,  which  correspond  to  air  and  argon,  respectively. 
The  numerical  results  have  been  correlated  to  within  3  percent  for  W  >  1.5 
and  to  within  2  percent  for  W  >  2.  0  by  the  expression 


F  = 


W(W  -  1) 

+  1. 25W  -  0.  80 


(17) 


Substitution  of  Eq.  (15)  into  Eq.  (5),  to  eliminate  i^,  yields  the 
following  value  for 


1.25W  -  0.80 
W(W  .  1) 


4/5 


(18) 


This  expression  has  been  used  with  the  values  of  in  Table  2  to  obtain 
corresponding  values  of  p^.  These  values  are  presented  graphically  in 
Fig.  4a. 

C.  NUMERICAL  RESULTS  FOR  i  ,  Re  ,AND  26  /d 

mm  m 

The  values  of  P^  in  Fig.  4a  have  been  used  to  compute  from  Eq.(4). 
The  results  are  given  in  Fig.  4b,  Standard  conditions  were  assumed  to  be 
T^^  =:  522°R  and  p^^  =  76  cm  Hg  so  that 

(£^)  =  6.93  X  10°  ft"^  air  (19a) 


=  7.  39  X  10^  argon  (19b) 

In  order  for  the  present  theory  to  be  valid,  it  is  necessary  that  the 
boundary  layer  be  turbulent.  This  will  be  the  case  if  the  Reynolds  number 
at  considerably  larger  than  the  transition  Reynolds  number.  An 
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Table  1.  Numerical  Constants  Used  to  Obtain  Results  in  Table  2  and  Pigure 


Table  2.  Values  of  (Eq.  7) 


Case  I 

Case  II 

_ 

Case  III 

Case  IV 

M 

P 

M 

P 

M 

p 

M 

p 

8 

8 

s 

o 

5 

1.41 

0.  0263 

1.  35 

0.  0253 

6 

0.  0283 

6 

0.  0270 

1.  73 

b.  0326 

1.  58 

0.  0332 

8 

0.  0220 

8 

0.  0210 

2.  24 

0.  0345 

1.  89 

0.  0389 

10 

0.  0188 

10 

0.  0175 

3.  00 

0.  0312 

2.  24 

0.  0413 

12 

0.  0157 

12 

0.  0149 

3.  46 

0.  0284 

2.  65 

0.  0413 

14 

0.  0129 

14 

0.  0122 

4.  12 

0.  0245 

3.  16 

0.  0391 

16 

0.  0116 

16 

0.  0106 

5.  20 

0.  0195 

3.  87 

0.  0348 

18 

0.  0104 

18 

0.  0095 

6.  08 

0.0164  ! 

5.  00 

0.  0282 

20 

0.  0094 

20 

0.  0084 

7.  55 

0.  0126  ! 

5.  70 

0.  0248 

-- 

— 

-- 

— 

8.  77 

0.  0105 

6.  71 

0.  0208 

-- 

— 

-- 

— 

10.  82 

0.  0080 

8,  37 

0.  0161 

-- 

— 

-- 

— 

— 

_ Z _ 1 

9.  75 

0.  0133 

-- 

— 

-- 

— 

(a)  Obtained  from  Eq.  (18)  and  Table  2 

Fig.  4.  Quantities  DeAning  Maximum  Separation  Distance  for  Turbulent  Wall 
Boundary  L^yer.  Cases  I  to  IV  Identified  in  Table  1 
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.18. 


appropriate  Reynolds  number,  based  on  i  ,  is  Re  =  u  (W  -  l)^i  /v 

m  me,o  me,o 

A  form  convenient  for  calculation  is 


This  parameter  is  plotted  in  Fig.  4c.  The  transition  Reynold's  number  has 
been  investigated  (Ref.  7)  and  is  in  the  range  0.  5  <  Re^  X  10”^  <  4  for 
1  <  Ms£9.  It  appears  to  increase  markedly  for  larger  M^,  due  to  the 
stabilizing  effect  of  the  low  wall  temperatures,  but  only  limited  dataA.re 
available.  The  present  turbulent  boundary  layer  theory  should  be  applicable 
if  Is  significantly  larger  than  Re^  (say,  >  5  Re^). 

It  has  also  been  assumed  (e.  g.  ,  Eqs.  1,  8)  that  the  boundary  layer 

thickness  is  small  relative  to  the  tube  radius.  This  condition  will  be  satisfied 

if  26  /d  is  small.  Here,  6  is  the  boundary  layer  thickness  at  i  .  If  6 
m  m  T  1  mm 

is  assumed  to  vary  as  in  the  case  of  a  uniform  external  stream,  it  can  be 
shown  that 


-  i  /I  -  w\ 

d  2  W  -  1  \6*/6  / 


(21) 


This  parameter  is  plotted  in  Fig.  4d.  It  is  seen  that  26^/d  >  1  for  <  2 
so  that  the  assumption  of  a  small  boundary  layer  thickness  is  violated.  For 
argon  considered  as  an  ideal  gas  (Case  I),  the  quantity  doesn't  go 

below  0.6;  and  for  ideal  air  (Case  II),  it  doesn't  go  below  0.45.  The  real  gas 
solutions  for  air  (Cases  III,  IV)  have  values  below  0.  5  for  >  6.  However, 
it  should  be  noted  that  26^/d  provides  a  very  conservative  estimate  for  the 
validity  of  the  present  theory.  What  is  mainly  required  is  that  the  boundary 
layer  not  become  fully  developed  (i.e.  ,  26^/d  <  1)  and  that  the  integrand  in 
Eq,  (8)  not  require  the  factor  'rr(d  -  2y)  (for  circular  cross  sections)*  The 
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second  of  these  conditions  is  probably  satisfied  when  26^/d  <  1  since  most  of 
the  boundary  layer  excess  mass  flow  occurs  in  the  portion  of  the  boundary 
layer  near  the  wall.  For  example,  by  neglecting  density  effects  and  assuming 
a  l/7th  power  velocity  profile,  it  can  be  shown  that  56  percent  of  the  excess 
mass  flow  occurs  within  a  distance  from  the  wall.  The  higher  gas 

density  at  the  wall  (particularly  for  strong  shocks)  will  considerably  increase 
the  fraction  of  excess  mass  flow  in  this  distance.  The  present  approach  is 
therefore  considered  probably  valid  for  >  3. 


III.  VARIATION  OF  SEPARATION  DISTANCE  WITH 
DISTANCE  FROM  DIAPHRAGM 


In  Section  II,  the  limiting  value  of  the  separation  distance  between  the 

shock  and  contact  surface  was  found  as  a  function  of  M  .  This  value  defines 

s 

the  test  time  in  shock  tubes  that  are  sufficiently  long  to  permit  the  separation 
limit  to  be  reached.  In  order  to  estimate  test  times  in  shorter  shock  tubes, 
it  is  necessary  to  consider  the  rate  at  which  separation  distance  increases 
with  distance  of  the  shock  from  the  diaphragm.  This  problem  is  now 
discussed,  using  the  approach  in  Ref.  4.  It  will  be  assumed  that  the  shock 
moves  with  uniform  velocity  and  that  the  flow  between  the  shock  and  contact 
surface  is  steady  in  a  shock- fixed  coordinate  system. 

A.  SEPARATION  BETWEEN  SHOCK  AND  CONTACT  SURFACE 

Consider  flow  in  a  coordinate  system  in  which  the  shock  is  fixed  and 
the  wall  moves  with  velocity  u  .  Assume  that  at  t  =  0  the  contact  surface 
coincides  with  the  shock,  and  at  some  later  time,  t,  the  portion  of  contact 
surface  that  is  in  the  free  stream  is  located  at  i  (Fig,  5).  The  problem  is  to 
find  i  as  a  function  of  t.  (In  the  present  section,!  is  considered  a  dependent 
variable  defining  the  location  of  the  contact  surface  relative  to  the  shock  as 
a  function  of  time.  ) 

For  steady  flow,  the  mass  flow  through  the  shock  equals  the  mass  flow 
through  a  control  surface  at  i.  This  yields  an  equation  identical  to  Eq.  (11), 
namely, 


“5  = 


P  u 
^e  e 

(p  u  ) 
e'o 


(11) 


where  T  is  again  the  ratio  of  the  excess  boundary  layer  mass  flow  at  i  to  the 
mass  flow  through  the  shock.  For  the  present  purposes  it  can  be  assumed 
that  the  excess  mass  flow  in  the  boundary  layer  varies  as  Then, 
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Fig.  5.  Flow  Model  for  Finding  Separation  Distance  Between 
Shock  and  Contact  Surface  as  a  Function  of  Time 


24. 


'  m' 


1-n 


Also,  it  will  be  assumed  that  p  /p  =1  (valid  for  strong 

e  e  I  o 


shock).  These  assumptions  have  been  shown  to  be  accurate  (Ref.  4),  except 
near  W 

(11)  then  becomes 


1  where  the  present  solution  for  is  invalid  anyway.  Equation 


(^) 


1  -n 


=  1  - 


u 


(22) 


e,  o 


Introduce  the  nondimensional  variables 


u  t  X 

r  e,0  S 

■  “jf  ■  wl 

m  m 


(23a) 


T  =  ^ 


m 


(23b) 


where  x  =  u  t  is  the  distance  of  the  shock  from  the  diaphragm  in  laboratory 

8  inf 

coordinates  (Fig.  6).  Noting  that  di/dt  =  u^,  write  Eq.  (22)  in  the  form 


T^*"=  1 


dT 

dX 


(24) 


Integrating, 


dT 


1  - 


T 


(25) 
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(a)  t-x  Diagram 


Fig.  6.  Flow  in  Laboratory  Coordinataa  Aaauming  Shock  Has 
Uniform  Velocity 


26 


Thus, 


-  ^  =  In  (1  -  t")  +  for  n  =  1/2  (26a) 

=  I  (in  i— I-  .  2  tan"^  t"  +  4T")  for  n  =  1/5  (26b) 

®  '  1  +  T*"  / 

Equation  (Z6a)  was  first  derived  by  Roshko. 

The  variation  of  T  with  X  is  tabulated  in  Table  3  and  is  plotted  in  Fig.  7 
for  n  =  1/2,  1/5,  Note  that  T  =  X  corresponds  to  an  ideal  (no  boundary 
layer)  shock  tube.  The  departure  of  the  curves  in  Fig.  7  from  T  =  X  repre¬ 
sents  the  departure  from  ideal  performance.  For  X  ^  0,  1,  the  departure  is 
small,  and  for  X^  10,  the  separation  distance  equals  the  limiting  value. 

Note  that  i  equals  59  percent  of  the  ideal  value  when  X  =  1,  n  =  1/5. 

B.  TEST  TIME 

In  the  previous  subsection,  the  separation  distance  was  obtained  as  a 
function  of  t.  A  quantity  that  is  perhaps  of  greater  interest  is  the  test  time 
(i.  e.  ,  the  difference  in  time  between  the  arrival  of  the  shock  and  the  arrival 
of  the  contact  surface)  at  a  fixed  distance  from  the  diaphragm,  x^.  This 
quantity  will  now  be  discussed. 

Designate  the  test  time  by  t.  For  x-^oo,  t-t  =i  /u.  Define 
®  ^  s  m  m  w 

T  -  t/t  ,  which  is  the  test  time  at  x  divided  by  the  test  time  at  x  oo. 
m  s'  s 

If  it  is  assumed  that  the  shock  moves  with  constant  velocity,  X  and  t  are 
related  (Ref,  4)  by 


X  =  X^  -  (T^/W) 


(27a) 


Equation  (26b)  was  derived  by  Dr.  E.  F.  Brocher  in  an  as  yet  unpublished 
study  of  test  time  for  the  case  of  small  departures  from  the  ideal  flow. 
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Table  3.  Variation  of  T  with  X  (£q.  25) 


T 

X 

n  =  1/2 

n  =  1/  5 

0.  01 

0.  0107 

0.  0101 

0.  05 

0.  0590 

0.  0527 

0.  10 

0.  128 

0.  110 

0.  20 

0.291 

0.238 

0.  30 

0.491 

0.388 

0.40 

0.  737 

0.  564 

0.  50 

1.  04 

0.  776 

0.  60 

1.43 

1.  04 

0.  70 

1.95 

1.39 

0.  80 

2.  71 

1.88 

0.90 

4.  04 

2.73 

0.  92 

4.  48 

3.  01 

0.  94 

5.  04 

3.  37 

0.  96 

5.  84 

3.87 

0.  98 

7.  22 

4.  73 

0.  99 

8.  60 

5.  60 

(J 


29 


Separation  Distance  as  a  Function  of  Distance  from  Diaphragm 
(TaUe  3) 


(27b) 


where  are  corresponding  values  from  Table  3.  The  variation  of  t 

with  X  is  plotted  in  Fig.  8a  for  W  =  2,  4,  8,  oo. 

An  analytical  expression  relating  T  and  f  can  be  found  as  follows. 

Assume  that  the  contact  surface  velocity  is  constant  during  the  period  of  time 

between  the  arrival  of  the  shock  and  the  arrival  of  the  contact  surface  at  x  . 

s 

The  test  time  would  then  equal  t  =  i/(u^  -  u^),  or,  from  Eq.  (22), 


WT 

T  =  - r— 

W  .  1  +  T 


(27) 


At  the  start  of  the  motion,  T  ~  WT/(W  -  1).  After  long  times,  t  =  T  =  1.  For 
large  W,  t  =  T  for  the  entire  motion. 

The  parameter  T  compares  the  test  time  t  to  the  asymptotic  test  time 
i  /u  .  It  is  also  of  interest  to  compare  t  to  the  test  time  for  an  ideal  shock 
tube,  T..  The  ideal  test  time  at  a  distance  x^  from  the  diaphragm  is 


T.  =  X  /(u  -  u  )W 

1  s  w  e.  o 


The  actual  time,  assuming  constant  between  arrival  of  shock  and  contact 

surface,  can  then  be  expressed  as 


JL  =  W  -  1 

^  W  -  1  + 


(28) 


This  is  plotted  in  Fig.  8b  for  n  =  1/5  and  W  =  2,  4,  8,  oo.  For  large  W, 
r/Xj,  =  T/X.  Note  that  when  X  =  1,  the  test  time  is  about  1/Z  of  the  ideal 
value.  The  ratio  t/t.  equals  (W  -  1)/WX  for  X  ^  5. 


-30- 


31 


-32 


(b)  Normalized  by  Ideal 


IV.  COMPARISON  WITH  ANDERSON'S  SOLUTION 


There  is  a  basic  difference  between  the  separation  distance  solution 

presented  by  Anderson  (Ref.  5)  and  the  one  presented  in  this  paper, 

Anderson's  method  for  finding  the  separation  distance  corresponding  to 

a  given  shock  location,  x^,  is  as  follows.  He  assumes  that  the  shock 

moves  with  uniform  velocity,  and  finds  the  separation  distance  by  equating 

Aip  ,  to  the  difference  between  the  ideal  amount  of  shocked  gas.  Ax  , 
*^e  o  8*^00 

and  the  amount  of  gas  that  has  slipped  past  the  contact  surface.  The  latter 

9/5 

he  finds,  incorrectly,  to  be  of  the  form  K^Ap^x^  (assuming  a  thin 
boundary  layer)  where  is  a  function  of  M^,  d,  p^,  and  is  tabulated  by 
Anderson.  The  separation  distance  is  then 


i  = 


t  is  seen  that  as  x  increases,  i  first  increases,  reaches  a  maximum  at  a 
finite  value  of  x  [namely,  x  =  (5/9  K|)^'  ],  and  then  decreases.  In  the 
present  solution,!  increases  monotonically  with  x^  and  asymptotically 
approaches  as  x^  -*’cx>*  Hence  the  two  solutions  have  different  behaviors. 

Anderson's  expression  for  the  net  mass  loss  past  the  contact  surface 
is  based  on  the  ideal  separation  distance  between  shock  and  contact  surface. 
It  can  be  readily  shown  that  his  estimate  is  too  large,  except  for  small  x^. 

If  the  free  stream  is  assumed  to  remain  uniform,  the  mass  flow  past  the 

4/5 

contact  surface  at  any  instant  is  proportional  to  i  .  The  net  mass  loss 
since  the  start  of  the  motion  for  a  shock  moving  with  uniform  velocity  is 
proportional  to 
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For  small  x  »  where  i  varies  linearly  with  x 

q/  c  ®  ® 

form  X  '  ,  as  derived  by  Anderson.  For  large  x  ,  i  approaches  i 

s  s  ..  m 

the  net  mass  loss  tends  to  become  proportional  to  i  '  x 

siderably  smaller  than  the  x^' ^  dependence  used  by  Anderson.  Thus, 

s 

Anderson's  solution  overestimates  the  mass  loss  and  underestimates  the 


the  net  mass  loss  has  the 

and 

which  is  con¬ 


test  time  for  large  x^.  This  accounts  for  the  differences  in  the  behavior  of 
the  two  solutions. 


5/4 

For  moderate  values  of  x^»  say  (5/9Kj)  ,  it  is  difficult  to 

predict  how  Anderson's  numerical  results  will  compare  with  ours.  Note 
that  the  present  theory  includes  the  interaction  between  the  boundary  layer 
and  the  inviscid  flow  solution),  which  reduces  the  predicted  test  time  as 
compared  with  a  non-interaction  theory  solution  and  Anderson's  solution). 
Anderson's  overestimation  of  mass  loss  tends  to  compensate  for  his  neglect 
of  the  boundary  layer-free  stream  interaction,  at  least  for  moderate  x^. 


In  Anderson’s  illustrative  examples,  M  =  6;  p  =10  cm  Hg; 

s  oo  ” 

d  =  3/2  in.;  x  =  15  ft;  and  d  =  4  in.;  x  =  50  ft.  His  theory  predicted  about 
s  s 

50  percent  of  ideal  test  time.  The  present  theory  predicts  about  30  percent 
of  ideal  test  time  for  these  cases.  Thus,  in  these  examples,  Anderson 
predicts  more  test  time  than  does  the  present  theory.  However,  with 
increase  in  x^,  his  theory  must  eventually  predict  lower  values. 
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V.  DISCUSSION 


The  present  results  can  be  used  to  estimate  shock  tube  test  time 
reduction  due  to  mass  loss  to  a  turbulent  wall  boundary  layer.  The  pro¬ 
cedure  is  as  follows*  First,  is  found  from  Fig.  4b  for  given  M^,  d, 

and  .  If  this  is  substantially  larger  than  Re^  (say  ^  5  Re^),  the 
limiting  separation  between  shock  and  contact  surface  will  be  due  mainly  to 
a  turbulent  boundary  layer.  Hence  the  turbulent  boundary  layer  analysis  can 

be  used  to  estimate  major  reductions  in  shock  tube  test  time.  Next,  i _  is 

m 

found  from  Fie.  4b  and  the  value  of  X  =  x  /Wi  can  be  determined  for  a 
^  s  m 

given  test  section  distance  from  the  diaphragm.  Finally,  the  test  time  can 
be  found  from  Figs.  8. 

If  R®^*  based  on  a  laminar  boundary  layer,  is  less  than  Re^>  laminar 
boundary  layer  theory  is  applicable  and  a  similar  procedure  is  followed 
using  the  corresponding  curves  in  Ref.  4. 

The  parameter  that  primarily  determines  R®^^  ^  given 

value  of  *  f^®  length-to-diameter  ratio,  x^/d,  determines  the  amount 
of  test  time  reduction.  Typical  values  for  these  quantities  can  be  obtained 
from  Figs.  9  and  10,  for  turbulent  and  laminar  flows,  respectively.  These 
Egures  are  based  on  the  identities 


(29a) 


(29b) 


where  the  terms  in  square  brackets  have  been  evaluated  for  both  laminar 
(Ref.  4)  and  turbulent  wall  boundary  layers  (Fig.  4). 
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Values  of  Re  are  plotted  versus  M  for  fixed  values  of  dp  in 
Figures  9a  and  10a;  the  product  dp^  has  the  units  (in.)  (cm  Hg).  Also 
plotted  in  these  figures  are  estimated  values  of  Re  and  5  Re  .  (It  was 
assumed  that  Re  varies  linearly  on  the  loglog  plot  from  Re  =  0.5  X  10^  at 

^  L  ^ 

Mg  =  1  to  Re^  =  4  X  10  at  =  9.  The  curves  are  not  continued  beyond 
Mg  =  9  because  of  uncertainty  in  Re^. )  Turbulent  boundary  layer  theory  is 
applicable  for  values  of  that  result  in  ^  5  These  values 

(from  Fig.  9a)  are 


dp 

dp 


air 

>  4  for  M  =3 

dp 

argon 

>  2  for  M 

s 

s 

>  10  for  M  =8 
s 

8 

cu 

>  10  for  M 

=  3 


=  8 

8 


Similarly,  laminar  boundary  layer  theory  is  applicable  when  Re^^  4  R®^» 
(from  Fig.  10a)  when 


air 

dp  <  0.3  for  M  =3 
*^00  s 

dp  <  0,  3  for  M  =8 
^oo  s 


argon 

dp  £0.6  for  M  =  3 
*^ao  s 

dp  <  1  for  M  =  8 


It  follows  that,  roughly,  turbulent  theory  applies  for  ^  3.nd  laminar 

theory  applies  for  dp^  £0.5,  for  3  <  Mg  £  9. 

In  Figures  9b  and  10b  Xg/Xd  is  plotted  versus  Mg  for  various 

The  ordinate  may  be  viewed  as  the  value  of  Xg/d  corresponding  to  X  =  1; 

that  is,  the  length->to->diameter  ratio  that  results  in  about  1/2  the  ideal 

test  time.  For  air  or  argon  turbulent  boundary  layers  and  Mg  >  3# 

X  /Xd»45,  80,  140  for  dp  5,  50,  500.  (Thus,  a  shock  tube  with 
s  oo 

dp^  »  5  and  Xg/d  »  50  will  have  about  half  the  ideal  test  time).  Similarly, 
for  air  or  argon  )aminar  boundary  layers.  Fig.  10b  shows  Xg/XdM  100,  10, 
and  1  for  dp^»0.5,  0.05,  and  0.005,  respectively. 
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Fig.  9.  Dependence  of  *nd  x^/dX  on  and  dp^  for 
(  )  Turbulent  Boundary  Layer 
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Fig.  10.  Dependence  of  Re  and  x  /dX  on  M  and  dp  for 
Laminar  Boundary  Layer  ^ 


-  AIR  (Pflo»0.076  CM  Hg,  REF.  4) 

- ARGON  (IDEAL  GAS,  REF.  4  ) 


2  4  6  8  10  12  14 


Ms 


(b)  xJdX 


Fig.  10.  (Continued) 


t 

These  numerical  results  are  intended  to  give  an  indication  of  when 
turbulent  or  laminar  theory  applies,  and  rough  estimates  of  the  test  time 
limitation.  More  accurate  estimates  are  made  from  Figs.  4  and  8, 

Before  the  present  generation  of  low  pressure  shock  tubes,*  the 
typical  low  pressure  section  of  a  shock  tube  had  an  internal  diameter  in  the 
range  1  in.  <  d  <  4  in.  and  a  length- to -diameter  ratio  in  the  range  40-150 
(Ref.  8).  Longer  tubes  were  sometimes  used,  particularly  in  studies  of 
shock  attenuation;  for  example,  a  2-in.  X  40-ft  tube  (x^/d  =  240)  was  used 
in  Ref.  9,  and  a  3. 75-in.  X  120-ft  tube  (x^/d  =  384)  in  Ref.  10. 

In  a  given  shock  tube,  the  largest  reduction  in  test  time,  due  to  a 

wholly  turbulent  boundary  layer,  occurs  at  dp^  »  5  (in.)  (cm  Hg),  for  which 

X  /Xd  »  45.  Tubes  with  x  /  d  «  45  will  have  X  =  1,  or  around  half  the  theo- 
s  s 

retical  test  time  (see  Fig.  8b).  Tubes  with  x^/d  «  150  will  have  X»  3,  or 

around  1/4  of  the  ideal  test  time.  (The  40  percent  lower  limit  on  test  time, 

noted  in  Ref.  6  and  discussed  in  Section  I,  is  therefore  not  sufficiently  low 

for  this  case).  With  regard  to  the  larger  values  of  x^/ d  considered  in 

Refs.  9  and  10,  it  might  be  noted  that  X  «  5. 5  for  x^/d  =  240  and  X»  8.  5 

for  X  /d  =  384.  Here,  the  test  times  are  approximately  0.  15  and  0.  1  of  the 
s 

ideal  values,  respectively.  For  these  values  of  X,  the  separation  distance 
has  reached  its  maximum  value  (Fig.  7). 

The  smallest  reduction  in  test  time,  due  to  a  wholly  laminar  boundary 
layer,  occurs  when  0.5.  From  the  results  given  in  Ref.  4,  it  can 

be  seen  that  the  test  time  is  about  1/2  the  ideal  value  when  x^/d  »  100.  The 
test  time  is  slightly  larger  than  the  corresponding  value  for  a  turbulent 
boundary  layer  with  dp  »  5. 


Present  day  low  pressure  tubes  have  low  pressure  sections:  6  in.  X  30  ft 
(AVCO),  24  in.  X  50  ft  (AVCO),  17  in.  X  36  ft  (Aerospace),  17  in.  X  70  ft 
(Caltech)* 
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It  is  thus  seen  that,  for  a  given  shock  tube,  operation  at  an  initial 
pressure  such  that  ^  3  results  in  a  local  minimum  in  the  maximum 
possible  test  time.  Higher  values  of  will  result  in  thinner  turbulent 
boundary  layers  and  therefore  more  test  time.  With  a  decrease  in  ^P^» 
the  onset  of  a  wholly  laminar  boundary  will  first  tend  to  increase  the  test 
time;  however,  with  further  decrease  in  ^P^>  the  boundary  layer  mass 
flow  will  increase  and  severe  test  time  limitations  will  occur. 
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VI.  CONCLUDING  REMARKS 


The  effect  of  wholly  laminar  and  wholly  turbulent  wall  boundary 

layers  on  shock  tube  test  time  has  been  discussed.  It  has  been  shown  that* 

for  given  shock  Mach  number,  the  maximum  test  time,  t  »  and  separation 

m 

distance,  i^,  depend  primarily  on 


"'^m  ~ ^  ^  ~  ^^Prr,  ^  ~ 

m  m  00  00 

,5/4  1/4  for  n  =  1/5 


Since  the  ideal  test  time,  t.,  is  proportional  to  distance  from  the  diaphragm, 

X  ,  it  also  follows  that 
8 


T 

m 


/T, 


i(dp  )''* 

X  '  *^oo 

s 


for  n  =  1/2 

for  n  =  1/5 


These  equations  show  that  the  maximum  test  time  is  proportional  to  p^  for 

laminar  flows  and  to  p^^^  for  turbulent  flows.  Hence  the  laminar  wall 

*^00 

boundary  layer  case  is  more  sensitive  to  p^  . 

Examples  of  previous  studies  have  been  noted  where  turbulent  wall 
boundary  layers  have  reduced  the  test  time  from  0.1  to  0.5  of  the  ideal 
value.  The  rule-of-thumb  (see  Section  1)  that  assumes  about  one  half  of  the 
ideal  test  time  to  be  actually  available  is  therefore  not  dependable.  The 
methods  outlined  herein  should  be  used  to  estimate  test  time  when  necessary. 
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It  should  be  kept  in  mind  that  the  present  calculations  actually  provide 
an  upper  limit  on  the  test  time.  It  is  also  necessary  to  consider  (a)  non¬ 
ideal  diaphragm  rupture  (which  promotes  interface  mixing  and  causes  non¬ 
ideal  flow  near  the  diaphragm)  and  (b)  flow  nonuniformity  between  the  shock 
and  contact  surface.  When  the  maximum  separation  is  reached,  the  flow 
nonuniformity  between  the  shock  and  contact  surface  can  be  estimated  from 
Eq.  (22)  by  assuming  steady  isentropic  flow  in  a  shock- fixed  coordinate 
system. 

It  should  also  be  noted  that  the  validity  of  the  turbulent  boundary  layer 
theory  used  in  the  present  report  has  not  been  established  for  strong  shocks 
(see  Appendix  A).  The  present  results  are  therefore  intended  as  a  first 
estimate  for  the  effect  of  turbulent  boundary  layers  on  shock  tube  test  time. 
Experimental  confirmation  of  these  results  would  be  desirable, 

t 
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APPENDIX  A 


Turbulent  Boundary  Layer  Behind  Moving  Shock 


( 


The  turbulent  boundary  layer  behind  a  moving  shock  was  analytically 
studied  in  Refs.  11  and  12.  The  general  validity  of  some  of  the  analytical 
results  given  in  Refs.  11  and  12  was  experimentally  established  in  Refs.  13 
and  7  for  Mach  numbers  up  to  about  3  and  6,  respectively.  Higher  Mach 
numbers  were  not  studied  therein;  hence  the  accuracy  of  the  results  in  Refs. 
11  and  12  has  not  been  fully  established,  particularly  for  strong  shocks 
(M^  >  6).  Nevertheless,  the  results  in  Refs.  11  and  12  for  displacement 
thickness  will  be  summarized  here  and  put  in  a  form  suitable  for  studying 
test  time  in  shock  tubes.  This  will  provide  at  least  a  first  estimate  of  the 
effect  of  turbulent  wall  boundary  layers  on  shock  tube  test  time. 

Assume  that  a  shock  moves  with  constant  velocity,  that  the  inviscid 
flow  behind  the  shock  is  uniform,  and  that  the  wall  boundary  layer  is 
turbulent  (Fig.  "2a).  The  variation  of  displacement  thickness  with  distance 
behind  the  shock  (from  Refs.  11,  12)  is 


(A-1) 


where 


(A.2) 


-45. 


Here,  6  is  the  boundary  layer  thickness,  and  9  the  momentum  thickness. 
The  subscript  m  refers  to  conditions  evaluated  at  a  mean  reference  static 
enthalpy  defined  by 


m 


e,  o 


{A.3) 


where  h^  is  the  recovery  enthalpy  (i.e.  ,  the  wall  value  of  h  for  which  there 
is  no  heat  transfer). 

In  order  to  find  K  ,  it  is  necessary  to  evaluate  {6>9‘/6)/{l  -  W)  and 
o 

(9/6)/ (1  -  W),  which  are  functions  of  W,  h  /h  ,and  h  /h  that  can 

r  e  I  o  Av  e ,  o 

be  found  from 


6»/S  _  1  j  ^e, o 

1  .  W  "  W  -  1  h^ 


[w  -  (W  -  l)t] 

- - — 

1  +  b;  - 


dt  - 1 


(A-4a) 


e/6 

1  -  W 


=  7 


w 


/ 


4  V 


a;  ,A.4k, 


1  +  b4  .  Ct* 


where 


(A-5) 


[in  reviewing  the  work  of  Refs.  11  and  12,  it  appears  that  for  strong  shocks 

it  might  be  advisable  to  replace  h  /h  and  h  /h  in  Eqs.  (A-4)  and  (A-5) 

e ,  o  V/  r  vr 

W  P  /p  P  /p  »  respectively.  However,  this  point  will  not  be 

w  c  s  o  r 

pursued.]  Equation  (A-4)  can  be  evaluated  for  given  b,  c,  and  W,  using  the 
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tables  presented  in  Refs.  10  and  11.  For  strong  shocks,  these  tables 
require  taking  the  differences  of  large,  nearly  equal  numbers;  thus  the 
accuracy  is  reduced.  (Approximate  expressions  for  evaluating  Eqs.  (A»4) 
for  strong  shocks  are  given  in  Appendix  B.)  If  convenient,  it  is  preferable 
to  evaluate  b^/6  and  0/|i  by  the  numerical  integration  of  Eq.  (4), 

Ideal  Gas 

Eqs.  (A«-4)  and  (A-2)  can  be  evaluated  as  a  function  of  for  gases 

with  constant  specific  heat  ratios,  y,  obeying  the  perfect  gas  laws,  p  =  pRT 

and  h  =  c  T.  For  these  gases, 

P 


e,o 


1  + 


(W  -  1)^ 
ZW  -  1 


r(0) 


(A-6a) 


e,  o 


W(Z  -  W) 
ZW  -  1 


(A-6b) 


1  /3 

where  Z  =  (y  +  l)/(y  -  1),  and  r(0)  =  (cr)  '  is  the  recovery  factor 

(or  =  Prandtl  number).  Eq.  (A-6a)  follows  from  the  definition  of  recovery 

factor,  and  Eq.  (A-6b)  is  a  normal  shock  relation.  The  latter  assumes  that 

the  wall  temperature  is  the  same  as  the  gas  temperature  upstream  of  the 

shock. 


Equations  (A-4)  have  been  evaluated  for  y  =  7/5,  r(0)  =  0.897  and  for 
y  =  5/3,  r(0)  =  0.875.  These  values  correspond  to  ideal  air  and  to  ideal 
argon,  respectively.  The  cases  y  =  7/5,  r(0)  =  1  and  y  =  5/3,  r{0)  =  1  have 
also  been  evaluated  to  find  the  effect  of  r(0).  The  results  are  given  in 
Fig.  11  and  can  be  approximated  by  the  following  expressions. 
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M-l  .  M-l 

«/e 


48- 


e/6  .  W  B  .  W  4  (7/3) 

1  -  W  *  C  “  80/3 


for  \  =  7/5 


(A.7a) 


W  -I-  2 
24 


for  \  =  5/3 


{A-7b) 


6#/6  ^ 
1  -  W 


D  =  0.  157 


for  y  =  7/5,  W  >  2  (A-7c) 


=  0. 176 


for  V  =  5/3,  W  >  2  (A-7d) 


These  expressions  agree  with  the  numerical  integrations  to  within  3  percent. 

Equations  (A -7)  have  been  used  to  evaluate  as  a  function  of 
using  the  Suthe^^  nd  viscosity  relation  for  argon  and  for  air.  The  results 
are  given  as  Cu-bci*  I  and  II  in  Table  2  in  terms  of 

Real  Gas 

For  strong  shocks  the  previously  noted  ideal  gas  relations  are  no  longer 
valid.  The  following  procedure  has  therefore  been  used  to  obtain  an  estimate 
of  the  turbulent  boundary  layer  behind  a  strong  shock  in  air  (M^  ^  6). 

The  problem  is  to  evaluate  (Eq.  A-2)  for  a  strong  shock.  First 
it  was  assumed  that  Eqs.  (A- 7 a)  and  (A -7c)  could  be  used  with  W  obtained 
from  the  equilibrium  shock  solutions  of  Ref.  14.  [Eqs.  (A-7a)  and  (A-7c) 
agree  with  the  limiting  values  of  6/6  and  6^/6  obtained  by  strong  shock 
assumptions*  Eqs.  (B-3),  to  within  5  percent  and  8  percent  for  W  =  10  and 
15*  respectively.]  It  then  only  remains  to  evaluate  the  term  [  ]  in 

Eq,  (A-2).  The  energy  equation  for  a  strong  shock  gives 


e*  o 


1  + 


W  -  1 
W  +  1 


r{0) 


(A.8) 


with  r(0)  =  0.897  for  air.  Reference  14  was  used  to  find  W  and 

X  /T^  =  T  /T^  for  given  M  ,  p^,  and  T  =  300®K.  It  was  then 
e*o  00  e,o  w  **  s  oo  w 
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assumed  that  p  /p  w  h  /h  «  T  /T  .  The  latter,  with  Eq.  (A-8), 

^e,o  w  e,o  w  e,o  n  \ 

permits  h  /h  »  T  /T  »  p  /p  to  be  found  and  Eq,  {A-2)  to  be 
me,o  me,o  e,o  m 

evaluated.  The  assumption  that  h  »T  1/p  is  only  approximately  correct 
but  should  be  sufficiently  accurate  for  the  present  purposes  (i.e.  ,  evaluation 
o£[  ]  in  Eq.  A-2). 

The  resulting  values  of  are  also  given  in  Table  2  (in  terms  of 
This  table,  together  with  Eq,  (7),  defines  the  turbulent  boundary  layer 
behind  a  strong  shock  in  air.  It  should  be  kept  in  mind  that  the  turbulent 
boundary  layer  theory  of  Refs.  11  and  12  has  not  been  verified  for  strong 
shocks  and  that  further  simplifying  assumptions  were  made  to  obtain  K^. 

Boundary  Layer  Theory  for  Local  Similarity  Approximation 

The  turbulent  boundary  layer  equations  will  be  put  in  a  form  that  is 
appropriate  for  an  origin  at  i. ,  a  wall  velocity  u^,  and  free  stream  conditions 
u^,  p^,  which  are  isentropically  related  to  conditions  directly  behind  the 
shock  (u  ,  p  ).  This  boundary  layer  is  illustrated  in  Fig.  3  and  is 

6^0  C I o 

used  in  the  local  similarity  approximation  in  the  body  of  the  report. 

For  the  boundary  layer  in  Fig.  3,  Eq.  (A-1)  becomes 


/  V 

-6#  =  k| - - — )  (i  - 

\“w  -  %  /  ^ 


=  K 


(A-9) 


where  V  =  u  /u  .  Equations  (A-7)  have  the  form 
e  e  I  o 

9/6  _  B  +  (W/V) 

1  .  (W/V)  "  C 


(A- 10a) 


1  -  (W/V) 


(A- 10b) 


-50. 


D 


where  B,  C  and  D  are  constants.  If  the  free  stream  Mach  number  relative 

to  the  shock,  M  =  (u  /a  )  ,  is  small  (i.e.  ,  strong  shock),  then  the 
e ,  o  e  e  o 

variation  of  fluid  state  properties  in  the  free  stream  can  be  neglected 

(p  =  p  ;  p  =  p  ),  and  Eqs.  (A-9)  and  (A- 10)  then  yield 
e  e,o  e  eiO 


6 


4V(.6*) 

d 


He(i 


4/5 


(A- 11a) 


where 


(A- 11b) 


Note  that  B  =  7/3  for  air  and  2  for  argon  (Eq.  A-7), 


f 
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II 


APPENDIX  B 

Limiting  Forms  for  0/6  and  6*/6 


Limiting  forms  for  0/6  and  6’!'/6  are  presented  herein  for  the  cases  of 
both  very  weak  shocks  and  very  strong  shocks. 

Weak  Shock 

For  the  case  of  a  weak  shock,  Eqs.  (A-4)  become 

* 

=  J[1  +  0(W  -  1)1  (B-la) 

=  1  +  0(W  -  1)]  (B-lb) 


Strong  Shock 

For  strong  shocks,  b»l,  c»  1  and  the  ratio  c/b  has  a  limiting  value 
given  by 


c  _  r(0) 

F"  7TTT5) 


Ideal  Gas  (B-2a) 


_ rJO) _ 

[(W  +  1)/(W  -  1)]  +  r(0) 


Real  Gas  (B-2b) 


Thus  c/b  =  0.  391,  0.  344  for  y  =  1-4,  r(0)  =  0.897  and  y  =  5/3,  r(0)  =  0.875, 
respectively.  These  are  the  values  for  ideal  air  and  argon.  For  a  real  gas 
with  r(0)  =:  0.9,  it  follows  that  c/b  =  0.  375,  0.424,  and  0.440  for  W  =  5,  10, 
and  15,  respectively. 


t 
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If 


terms  of  order  b 


-1 


compared  with  1  arc  neglected »  Eqs.  (A-4) 


yield 


_  1 
’  W  -  1 


Tw[rTT)  *  ']l  ‘ - 

N=0 


{B-3a) 


e/5 
1  -  w 


+  0(b’Sl 


(B-3b) 


The  constants,  B,  C,  D,  E  have  been  evaluated  for  several  values  of  c/b 
and  are  given  in  Table  4. 

For  strong  shocks  in  air,  Eqs.  (A-Ta)  and  (A-7c)  agree  vidth  Eqs. 
(B-3)  to  within  5  percent  for  W  =  10  and  to  within  8  percent  for  W  =  15. 
Hence  Eqs.  (A- 7a)  and  (A- 7c)  are  sufficiently  accurate  for  a  first  estimate 
of  0  and  behind  strong  shocks  in  air. 

Table  4.  Constants  Defining  9/6  and  6^/6 
for  Strong  Shocks  (Eqs.  B*3) 


c/b 

B 

C 

D 

D 

WSM 

3.  16 

28.2 

WBM 

BB 

mssM 

3.  19 

29.  2 

mmm 

0.  45 

3.23 

30.  3 

mmm 

Hg 

0.  50 

3.27 

31.6 

HH 

Hi 
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